UAB-FT-643 



DiRAC Vs. Majorana Neutrino Masses 
From a TeV Interval 



D. Diego 1 and M. Quiros ^ " 

" Theoretical Physics Group, IFAE/UAB 
E-08193 Bellaterra (Barcelona) Spain 

Institucio Catalana de Recerca i Estudis Avangats (ICREA) 



Abstract 

We investigate the nature (Dirac vs. Majorana) and size of left-handed neutrino masses in a 
supersymmetric five-dimensional model compactified in the interval [0,7ri?], where quarks and 
leptons are localized on the boundaries while the gauge and Higgs sectors propagate in the bulk 
of the fifth dimension. Supersymmetry is broken by Scherk-Schwarz boundary conditions and 
electroweak breaking proceeds through radiative corrections. Right-handed neutrinos propa- 
gate in the bulk and have a general five- dimensional mass M, which localizes the zero modes 
towards one of the boundaries, and arbitrary boundary terms. We have found that for generic 
boundary terms left-handed neutrinos have Majorana masses. However for specific boundary 
configurations left-handed neutrinos are Dirac fermions as the theory possesses a conserved 
global f/(l) symmetry which prevents violation of lepton number. The size of neutrino masses 
depends on the localization of the zero- modes of right-handed neutrinos and/or the size of the 
five- dimensional neutrino Yukawa couplings. Left-handed neutrinos in the sub-eV range re- 
quire either MR ~ 10 or Yukawa couplings ~ 10^'^ i?, which make the five-dimensional theory 
perturbative up to its natural cutoff. 
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1 Introduction 



Within the Standard Model (SM) framework the origin, nature and hghtness of neutrinos are 
still open questions. A commonly accepted explanation for Majorana neutrinos is the so called 
seesaw mechanism [1, 2]. It requires the existence of a sterile ^ and massive (M^ ^ Mz) 
right-handed (RH) spinor weakly coupled to the SM left-handed (LH) neutrinos and the Higgs 
through a Yukawa term yielding thus a light Dirac mass {mo -C Mj{). The lowest mass 
eigenvalue is then mj, ~ -g^ and it is ~ meV if the Yukawa coupling is weak enough and/or 
the RH neutrino mass is large enough. It was originally proposed in the context of four- 
dimensional (4D) models and it predicts the existence of new physics at the scale Mr. The 
drawback of the seesaw mechanism is that it either requires an extremely large value of [out 
of reach for the Large Hadron Collider (LHC)] or introducing a small value of the 4D Yukawa 
coupling ~ 10^^. The alternative possibility that neutrinos are Dirac fermions would require 
much smaller 4D Yukawa couplings ~ 10~^^. In any way the nature of neutrinos (Dirac or 
Majorana), which is essential for many experimental signatures, still remains as the big open 
question in neutrino physics. 

Some of the above problems may find a solution in the context of theories with (compact- 
ified) extra dimensions. They provide a natural scale, the inverse compactification radius of 
extra dimensions 1/R, and moreover they can "naturally" provide very small 4D Yukawa cou- 
plings if some of the fields participating in the localized Yukawa interactions are exponentially 
localized far away. This led to the physical scenario where the RH neutrino belongs to a hy- 
permultiplet propagating in the bulk of the extra dimensions while Yukawa interactions are 
localized at fixed points of the bulk, e.g. in orbifold compactifications. 

In general in a theory defined in an extra dimensional scenario two main points should be 
answered: 

• Why neutrinos are so light? 

• What is the nature of neutrinos: Dirac vs. Majorana fermions? 

Although lots of different studies have been done so far [3]- [23] in supersymmetric or non- 
supersymmetric theories, and flat or warped space-time a systematic analysis of the previous 
^With respect to the SM gauge group. 
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questions was not addressed in a compactified theory with arbitrary boundary terms. 

In this paper we have investigated in detail the hghtness and nature of neutrinos in a 
supersymmetric five-dimensional (5D) theory defined on the interval x [0, irR] where gauge 
and Higgs bosons propagate in the bulk of the fifth dimension and SM fields belong to chiral 
superfields localized in one of the boundaries of the 5D space. We will also consider the 
compactification scale in the TeV range, which can lead to production of KK excitations [24, 25] 
at the LHC ^ and where supersymmetry can be globally broken by the Scherk-Schwarz (SS) 
mechanism [26]. In these theories the supersymmetric spectrum has a very characteristic 
pattern [27], which can very easily be identified in experiments, and electroweak symmetry 
breaking (EWSB) can proceed by radiative corrections [29, 30] bypassing some previously 
outlined difficulties in Ref. [28]. 

The plan of the paper is as follows. In section 2 the most general supersymmetric action for 
RH neutrinos in the bulk is presented as well as the boundary terms, and boundary conditions 
arising from the variational principle are written. In section 3 the wave functions for the 
bosonic and fermionic components of the RH neutrino hypermultiplet satisfying the equations 
of motion are deduced and the number of degrees of freedom identified. In section 4 the 4D 
theory obtained after integrating the fifth dimension is worked out and the effective action for 
the light neutrino is computed in the background of the Higgs field. Our main physical results 
can be found in section 5 where the different situations concerning the nature and size of LH 
neutrino masses are discussed. Finally section 6 contains our conclusions. 

A short summary of the main results is appended now. Concerning the nature of the 
LH neutrino mass we have found that everything depends on the relative orientation of the 
boundary action (characterized by a vector s) with respect to the bulk action of the RH 
neutrino characterized by a mass vector p. In the generic situation where both vectors s 
and p are arbitrarily oriented, lepton number is violated and we find a Majorana mass for 
LH neutrinos. In the particular case where s and p are parallel (or anti-parallel) there is a 
conserved global U{1) symmetry group which prevents lepton number breaking and produces 
Dirac masses for neutrinos. As for the size of the LH neutrino masses, in the case of Dirac 
masses they can be in the meV range if the zero-mode of the RH neutrino is exponentially 
localized towards the brane opposite to that where the SM fermions localize. In particular if 
^This possibility will of course depend on the detailed localization properties of quarks and leptons. 
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M is the localizing mass of RH neutrinos the correct size for the LH neutrino mass requires 
that MR ~ 10. In the generic case where LH neutrino masses are Majorana, in order to get 
the correct size without imposing anomalously small 5D (dimensional) Yukawa couplings it is 
shown that quarks and leptons should be localized in opposite boundaries while the bulk Higgs 
localizes towards the quark boundary. In this case the Yukawa couplings should be somewhat 
small (~ 10~^R) although for particular values of s ■ p there is a cancellation in the neutrino 
mass matrix for large values of M and the correct spectrum of neutrinos can be again obtained 
for MR ~ 10. 

2 The action and the spectrum 

The action for the sterile neutrino is that of a 5D = 1 supersymmetric theory with boundary 
terms defined in the manifold S = M4 x J, where M4 is the Minkowski space and / is the interval 
[0,7ri?], R being the compactification radius that we assume to be ^ TeV^^. Thus the right 
handed neutrino belongs to a (SM singlet) hypermultiplet (A^, iVc) with A^(c) = 0(c) + \/29tJj(^c) + 



The explicit form of the action is 



where 



S — Shk + + S'bd , (2.1) 

^bk = / d^O [NN + N^N,]- I (fe Ncd^N + h.c., (2.2) 
Jt, Jt. 

h.c. , (2.3) 



T 



aN,N + -N^--N^ 
2 2 " 



5bd = J^^ d'e l^^N' + ^N^ + uNN^j + h.c. , (2.4) 

and /i, A and u are arbitrary (dimensionless) complex numbers ^, while the bulk mass terms 
only depend on two (mass) parameters (a G M, 6 G C) in order to guarantee the SU{2)r 
invariance. In fact the fermionic component of Eq. (2.3) provides the most general fermion 
mass Lagrangian invariant under the 5D Lorentz transformations 

a'ipci^ + ~ y^c^c + h.c = \1/ + vE''^ + — (2.5) 



^It is understood that boundary terms occur independently at y ~ and y = ttR. 
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^ = ( 1" ) (2.6) 



where \l/ is the Dirac spinor 

= Cs^"^ its charge conjugate and C5 = 727075 the 5D charge conjugation operator ^. 
By applying the variation principle on (2.2)-(2.4) we find the boundary conditions 

fxN + i^N^- N^ = 0, \N^ + uN = 0, (2.7) 

which turn out to be a set of manifestly supersymmetric boundary conditions. One can easily 
check that the system of equations (2.7) is overdetermined unless 

/iA-(z/-l)i/ = 0, (2.8) 

and that it is invariant under the redefinitions 

^ N, \ ^ u ^ 1 - u . (2.9) 

In the special case u = the boundary conditions read 

A = 0, fxN -Nc = 

or . (2.10) 

/i = 0, iV, = 

while the case = 1 is obtained from the previous one by means of the relations (2.9). In the 
general case where 7^ 0, 1 Eq. (2.7) reduces to 

^N^ + N = 0, (2.11) 

which means that the complex parameters {u, fi, A) are highly redundant since only the complex 
number z = X/u does matter. Notice that by letting z to take any complex value we cover the 
whole set of boundary conditions (including u = 0, which corresponds to ^ — 00). Actually 
one can easily show that the whole complex plane is covered consistently with (2.8) by the 
restriction 

ueR, fi = -X*. 
As a matter of fact a possible parametrization is given by 

^ = ^(1 + 53), /i=-A* = is_, (2.12) 



^Note the different definitions of the 5D (C5) and 4D (C4 = — ^7270) charge conjugation operators which 
gives rise to the minus sign which appears in Eq. (2.3). 
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with S- = Si — is2 where st are three real parameters such that s ^ = 1. 
The boundary conditions now read as 



N 



0. 



(2.13) 



where a are the Pauh matrices. In components Eq. (2.13) reads as 



2(1----) 



[l + s-a) [dy+p-sM] 



-(1-s-a) 
2 ^ ' 



respectively, where we have defined the shorthands 

1 









(2.14) 
(2.15) 



P 



ibR,-bj,a), M = Ja^ + \b\ 



The spectrum allowed by these boundary conditions can be read off from Refs. 
and it is provided be the zeroes of the function 



sin^(7rr) — (cq — Cn) —tan{'Kil R) 



cos {7TT ) + CoC-n 



tan2(7rfii?) 



with cos(27rr) = Sq ■ St^, Cf = p ■ Sf and = — M^, m being the physical mass eij 
In particular for cq = = c {r = 0) the spectrum predicted by (2.17) is 



mn = - VM^R^ + V? , n=l,2,3,--- 
R 



(2.16) 

[29, 30], 

(2.17) 
renvalue. 

(2.18) 



with = \ — (? 



3 Wave functions 

In this section we will explicitly write the solutions to the equations of motion and the boundary 
conditions for the model presented in the previous section. Let us first re-express the action 
in a more appropriate way. As one can see from (2.3) the most general mass term involves 
both Dirac (a) and Majorana (6, 6*) masses and there is a family of continuous transformations 
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parameterizing the possible mass configurations. These are SU (2) rotations acting in the space 
of chiral supermultiplets [NcN)"^ ^. Although those transformations do not leave the action 
invariant they are symmetries of the spectrum equation (2.17) ^. According to (2.16) the action 
(2.1) can be written as 

Ar^e(l-s-a)Ar| +h.c. , (3.1) 



1 

-'as 



where 

and 6 is the totally antisymmetric 2-tensor defined as 



^=('lr]^ (3-2) 



J)=i^2. (3.3) 

For simplicity we will be interested in a pure Dirac bulk mass term, i.e. p = (0, 0, 1). The set 
of unitary rotations bringing p to (0, 0, 1) is given by U{p) where 



with Cj; = p3, Sp- = y 1 - (ps) and 

e'^^ = (2sp^)"^/2 (^^/^Hri;^-i^^-p, + sp^ 

Notice that U{pi,p2, 0) is the limit of U{p) when p^ — ^ 0, i.e. 

Under these rotations the kinetic term remains invariant while the boundary matrices, like the 
bulk mass, transform covariantly. 



^Notice that this symmetry makes sense since superfields and Nc arc uncharged. 

^Actually the spectrum equation docs only depend on S'[/(2)-invariants. Geometrically speaking these 
(global) unitary transformations change the basis where the bulk and boundary matrices are expressed but the 
relations between them remain unaltered. 
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In the basis where the bulk mass is Dirac-hke we have sq = s^, 
of motion now read 



[—s,0,c). The equations 



((9^ - □ - M2) = 
{dl - □ - M2) 0, = 



for bosons, 



i a^d.i: - d,i.^_ - M^_. = 01 

ia^'^^^|;c + ^5^^J-M^P = j ^^rmions. 

3.1 BOSONIC SOLUTIONS 

For a given mass m the bosonic fields satisfy the Klein-Gordon equation □ 0(c) 
Therefore the general solution to the 5D equations of motion is given by 

$(a;, y) = A{x) cos{Qy) + B{x) sm{Qy) , 



(3.6) 



(3.7) 



-?^V(c)- 



(3.8) 



where $ = ( '^^ ] and f2 has been defined in (2.17). The boundary conditions (2.14) at y = 
impose the restrictions 

' 1 + c- s 



A 
B 



1 — c — s 
cM f l + c- s 



a{x) 



n \ i-c- s 



+ ( l + c-s ) ^^^^ 



where a{x) and b{x) are independent complex functions verifying the above Klein-Gordon 
equation. Finally the boundary conditions a.t y = ttR impose 



b{x) sin(fi7ri?) = , 



a(xi 



n 



sin(n7ri?) = , 



(3.9) 
(3.10) 



with two possible solutions 

1. sin(fi7ri?) 7^ and hence b{x) = and fi^ = — c^M^, whose eigenfunction is 



l + c-s 
1 — c — s 

l+c-s 
1 — c — s 



cM 

cos{ilRy) —sm{Q,Ry) 



n 



ip{x) 



e 



-McRy 



if{x) 



(3.11) 
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2. ilR = n G with eigenfunction 

+ 



1 + c- s 
1 — c — s 

-1+c+s 
1 + c- s 



(3.12) 



with r{y) = cosily) - ^sm{^y) and ^"(y) = sm{^y). 

3.2 Fermionic solutions 

Eq. (3.7) can be re-casted into a single Dirac equation as 

{irdf,-^'d,-M) * = 0, (3.13) 
with \l/ the Dirac spinor defined in Eq. (2.6), and whose formal solution is given by 

e(x), 



[^^^Q —M 

cos(V-n - IVP y) + 7^ / ^ = sin(V-n - IvP y) 



V-n - M2 

where 0(x) is the initial value at y = which fulfills the 4D Klein-Gordon equation □ B + 
m^O = 0. Thus the eigenfunction corresponding to the m-th mode reads 



cos(fiy) + 7 ^ - 



sin(f2i/) 



e(x). 



(3.14) 



Its initial value satisfies the boundary condition (2.15) at y = 0, that is: 

e = 



(3.15) 



(1 +c- s)x 
(1 -c- s)x 

with X ciii arbitrary 4D Weyl spinor. Since we are dealing with a 4D Dirac spinor, if ^ is the 
Dirac partner of x ^ satisfying the equations of motion 

ia''df,^=mx, ia^df,x = m^, (3-16) 
by plugging now (3.16) and (3.15) in (3.14) we obtain the fermionic wave eigenfunction 

(l+c-s)/f^(y)x 



+ m 



-il + c-s)^ 



sm{Qy) 



(3.17) 



^This is in straight analogy to the orbifold case where we have a Dirac spinor in the bulk although parity 
assignment projects out one of the components at the fixed points such that there we have a single Weyl 
(Majorana) spinor. 
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where 



f2{y)= cos{ny)±^sminy). (3.18) 



Finally the boundary condition at y = ttR imposes the further condition 
sm{il7rR) f 1 — c s \ f —M{l + c — s) m(l — c — s) \ f X 



Q \ s 1+c J \ M{1- c- s) -m(l + c-s) ^ V ^ ' ^ 
which is equivalent to 

0. (3.19) 



sin(n7ri?) f s M m \ fx 



Q \ sM m 
Eq. (3.19) again has two possible solutions: 

1. QR ^ Z with the solution 

m^ + sMx = 0. (3.20) 

The condition (3.16) implies now = X = V which corresponds to a Majorana spinor 
and 

m = -sM. (3.21) 

2. VLR = n with n G Z_|_. In this case we have two independent spinorial degrees of freedom, 
^" and x"", degenerated in mass. 

Thus we re-encounter the spectrum (2.18) and the corresponding wave functions turn out 
to be 



{l + c-s)f%y) r 
{l-c-s)fl{y)r 



Some comments about these solutions are now in order. On the one hand notice that 
they are not, in general, factorizable as f{y)g{x) which makes it explicit that the orbifold-like 
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ansatz is not suitable ®. On the other hand we want to remark that the presence of ^ does 
not appear as a unitarity problem. From (3.16) it can be expressed in terms of a^d^x which 
is in agreement with the uniqueness of the solution to (3.13). In fact ^ and x can be thought 
of as off-shell independent degrees of freedom. Finally from (3.14) and (3.6) one immediately 
checks that given a value of Q the solutions corresponding to ±Q are the same since there 
is no associated degeneracy. The sign of the root is a matter of convention since it can be 
absorbed by the redefinition y nR — y, which shows that both signs do correspond to the 
same eigenstate. 

From now on we will concentrate on the fermionic sector describing RH neutrinos. In 
particular we will write down the effective 4D action for the hyperfermions [ip, ipc) coupled 
to the SM sector which is localized at one of the boundaries, i.e. y = (/ = 0, tt). The 
Yukawa couplings will induce a Dirac mass connecting the left- and the right-handed neutrinos 
and we will find the eigenvalues (Dirac or Majorana) of the effective mass matrix in an extra 
dimensional generalization of the 4D see-saw mechanism ^. In addition we will discuss in detail 
the possibility of getting an ultra light mass eigenvalue with natural values of the 5D Yukawa 
couplings. 

4 Effective action for neutrinos 

In this section we will develop the 4D effective action for RH neutrinos coupled to the SM 
matter localized at the brane y = yj^. We will obtain the lowest mass eigenvalue in the neutrino 
sector by solving the characteristic polynomial of the infinite dimensional effective mass matrix 
involving the LH neutrino and the zero and non-zero modes of RH neutrinos in the background 
Higgs field. As we will show this can alternatively be done through the effective mass matrix 
for the LH neutrino and the zero mode of the RH neutrino resulting after integrating out the 
higher KK modes of RH neutrinos. For the Higgs field which propagates in the bulk and whose 
zero mode is exponentially localized towards one of the boundaries we will consider the action 
computed in Ref [30]. 

*In fact only in the particular case where s = 0, c = ±1 the solutions (3.22) break up into two orthogonal 
and factorizable functions. Furthermore in the M ^ Q limit both functions have a definite parity since the 
basis of mass eigenfunctions is in that case {cos(ny/i?), sin(nj//i?)}. 

®See Ref. [3] for a previous non-supcrsymmetric analysis. 
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The 5D action under study is 
1 



'eff 



2jE 
1 

2 

1 



(-^sV- + M^) + 7^ (95^/;, + Mij,) + h.C. 

^ + (S3 - l)V^c] + [(1 + Ss)^ - S+^e] + h.C. 

iPc'^lHc + h.C. 



(4.1) 



where denotes the LH neutrino, Y^, is the 5D Yukawa couphng (with mass dimension -1), 
t/fg stands for the boundary where is locahzed, i.e. /o = or /o = vr and is the lowest 
mode of the Higgs which according to Ref. [30] is given by He — ^/2Mh e~^'"y h{x) We 
now expand ip and ipc as 



'1 + c - s)e 



-Mcy 



+ 



. I (1 + c - s)r^{y) x"(x) + ^(1 - c - s) sin(ny/i?) ^(x) 

n>l 



^ ={l- c - s)e-^"^yr]{x) 



Y.\il-c- s)fl{y) x"(x) - 1^(1 + c - s) sHny/R) i 
L ^ n 



n>l 



and whose components satisfy the (free) equations of motion 



1 

1 



r 



(4.2) 
(4.3) 



with mn = VWWTn?. 

"'^"The Yukawa interaction terms are strictly localized on the boundaries since they are not SU{2)fi invariant. 

^^The theory of the Higgs hyperscalar {H, He) was worked out in detail in Ref. [30]. The bulk mass term 
in the Higgs sector Mh is further restricted by the SU{2) (g) U{1) invariance of the SM and there is a similar 
parameter to co,7r defined in section 2, ch, which is restricted by the condition of the Higgsino mass /x to be 
sh = h/Mh so that \ch\ — 1- 
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By integrating over the fifth coordinate one obtains the effective 4D action 



S, 



eflf 



n>l 

+ j d^xn V2MHe/, (1 + c-s) r/z/i/i 



where fco is the real number 



2cMttR 



/Cr, the tower of matrices 



fc„ = 27ri?^(l-. 



Mi? 

1 + S (Ti 



fn stands for 



and 



fn 



Xn 
in 



Q-MH\yf^-yfo\_ 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



where fh = OoTfh = TT depending on the boundary the zero mode is localized towards. By 
redefining the modes as 

1 



fco 



with 



l±s- 



MR 



nir, 



(4.8) 
(4.9) 

(4.10) 



12 



the 4D effective action (4.4) can be rewritten as 

1 



cS 



1 - 

2 



J2 d^x- {tp'^a^d^tP'^ + tpla^'df.^Pl) 

n>l 



n>l 



X 



+ h.c. 



where 



and 



y(o) = n(l + c-,),/^.;,.e,„ 



1 - S 



Yt^ = ^YAl+c-s)Mntf,kt^ 



are the 4D effective Yulcawa couphng constants. Here we have defined C/,^ as 



(4.11) 

(4.12) 
(4.13) 

(4.14) 



where fy depends on where the lowest mode of the RH neutrino locahzes towards. Once the 
Higgs gets its vacuum expectation value, {h) = v, the Yukawa couplings turn into Dirac mass 
terms and we are thus left with an effective mass matrix connecting LH and RH neutrinos 
given by 



1 



n>l 



2R 



zrrir 



(4.15) 



n>l 



Notice that higher modes of RH neutrinos appear in (4.11) as Majorana spinors albeit we 
started with Dirac fermions. However if we redefine the fields for n 7^ as 



(4.16) 
(4.17) 



we recover the expected Dirac spinors due to the degeneracy in mass of the higher modes. 
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Now we can find the eigenvalues of tlie infinite mass matrix (4.15) by computing its char- 
acteristic polynomial [3] 



/ -A 



P(A) = det 



m 



(0) 
D 



m 



mg^ sM-\ 



in) 
D+ 



m 



in) 
D- 



m 



m 



in) 
D+ 
in) 
D- 



—rrir, — A 



m„ — A 



V 



(4.18) 



where mg"* = vY'^^^ and m^^}^ = vY^\ The determinant of (4.18) yields 



P(A) 



(4.19) 



k>l 



A + 



m 



sM-X 



E 

n>l 



m 



in) 
D+ 



m 



in) 



m„ + A m„ — A 



When s 7^ the smallest eigenvalue, X^, will not be that canceling either A — s M or X^—m\. 
Therefore it should verify the equation 



Xl + 



m 



sM-Xl 



n>l 



m 



in) 
D- 



m 



in) 
D+ 



21 



TUn - Xl mn + Xl 



0. 



(4.20) 



Since m^^l_ and m^^' are ~ vY^/ R, and we are assuming this parameter to be much smaller than 



,(0) 



M, we can expand the solution for the dimensionless parameter Xl/M in powers of /5 — j^^j^ 



as 



Xl 
M 



(4.21) 



which makes sense whenever the lowest order is small. Substituting back in (4.20) we find that 
at lowest order A^ is given by 

2 



Xl 
M 



m 



R 



sM'^ ^ Mm. 

n>l 

2(l + c)s {tf,f RMhv'Y,^ 



m 



in) 
D~ 



m 



in) 
D+ 



TT 



E 



n 



n> 



^ (n2 + M^R'^){n^ + c^M^R^) 



(4.22) 
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Finally the series in (4.22) can be computed analytically by means of a Poisson re-summation 
giving 



A, 



-v^Y^^-^Mh {tf^f [2c(e/J^ + coth(7rMi?) - c coth(7rc Mi?)] . 



(4.23) 



Notice that the series in (4.22) is finite because both converge to the same limit when 
n oo. This feature is just a consequence of the structure of the bulk mass matrix in (2.3) 
which is consistent with 5D supersymmetry and Lorentz invariance. 

Alternatively this result can be obtained as the diagonalization of the effective action in- 
duced after integrating out the higher KK-modes ip"^ in (4.15) for momenta much smaller than 
their masses (i.e. neglecting their kinetic terms). From (4.15) the equations of motion for if)^ 
are given by 



Rm 



(n) 
D± 



rrir. 



(4.24) 



and substituting back in (4.15) we find, in matrix form, the following effective mass coupling 
for and ( = Uji 



(0) 

/i ml)' 
mg^ s M 



with 



E 

n>l 



R 



mr. 



m 



(n) 



(n) 

D- 



(4.25) 



(4.26) 



Since and are much smaller than M the light and heavy eigenvalues Xl,h of the mass 
matrix in (4.25) are 



A, 



m 



sM 



A 



H 



sM . 



(4.27) 



where the light eigenvalue coincides with that found in Eq. (4.23). 

A particularly interesting case is found when s = 0, that is when the boundary matrices 
are precisely aligned with the bulk mass matrix. In that case the characteristic polynomial 
simplifies to 



fo(A) 



mi- 



.k>l 



I 



n>l 



m 



D 



mi 



— m 



4°^ 



(4.28) 
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where 

^^S'f-2^HPM,H-^^^^^, (4,29) 

Notice that (4.28) is an equation for which means that both ±A are solutions and thus the 
set of eigenstates of the whole mass matrix are exactly degenerate by pairs and therefore they 
can be gathered to yield Dirac fermions. Following the same methods used above we find that 
the lowest eigenvalue is given by 

Xl± = ±2Y,v tf,^^/MHMe-^''yf-^ , (4.30) 

The effective mass matrix for pl, will be given in that case by 

7 V (4.31) 

where we have defined rriD = 2YyV e \/ MM^ C'^^^y^^ . 

The degeneracy of the spectrum for s = can be understood in terms of a symmetry which 
takes place only within this case. As a matter of fact s = means that the vectors p, Sq, are 
all aligned along the same direction and hence a f/(l) subgroup of unitary rotations around this 
direction leaves the action invariant. In terms of the fermion components these transformations 
translate into 

(77, x") - e^" (r/, x") , e-e-"r, z/L-e-^z/i (4.32) 

where a is a real parameter. Notice that this symmetry forbids any Majorana mass term 
in particular for i/^, and hence must vanish, as it is evident from Eq. (4.4). 

5 Discussion on neutrino masses 

In this section we will apply the previous results to discuss the possibility of getting, within 
this kind of models, an (ultralight) neutrino mass in the sub meV range. The first task will 
be to set the range of dimensional Yukawa couplings which appear in the 5D action in the 
leptonic sector 

^ {Y,,H,ULi^, + Y,HheB) (5.1) 



/ 



^This symmetry plays the role of the lepton number symmetry of the SM. 
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i.e. Yy^, with mass dimension -1, and Y^, with mass dimension -1/2, where I = {r, /i, e}. A 
naive estimate of wave function renormahzation corrections to the Yukawa couphngs in the 5D 
theory sets bounds as 

yi^^< -y= (5.2) 



R AR' - y/R ' VAi? 
so that taking AR ~ 20 we obtain 0{1) upper bounds on the dimensionless Yukawa couphngs 
Vuij- We can now distinguish two different scenarios where neutrino masses are either Dirac 
or Major ana: 



5.1 Dirac mass 

We will assume here that all the SM matter is strictly localized on the y = brane and the 
zero mode of the Higgs is localized towards it as well, thus ej^ = 1, while the zero mode of RH 
neutrino is exponentially localized towards y = nR, i.e. yf^ = tiR as Fig. 1 shows. 



Qli Ur, Dr 
LliEr 




y = ttR 



Figure 1: Bulk and brane matter distribution for a neutrino Dirac-like mass. 



We will consider for the moment the particular case where s = 0, i.e. where vectors p 
and so^,r are all aligned along the same direction. As it was shown in the previous section we 
obtain a Dirac mass connecting and z/r, which is of order ~ 2v Y^,^ y/ M Mh , with 
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e/j = 6 In Fig. 2 we show the Dirac mass as a function of MR for = 5 TeV, y^,^ ~ 1 

and MhR ~ 1.6 [30] We can see from Fig. 2 that m^^ < 1 eV for MR > 9 although m^^ 
decreases exponentially when MR increases and thus ^ 1 meV for MR ^11. In this 
scenario there is no wave function suppression for the charged leptons whose Yukawa couplings 
should therefore be given by ~ m^/v. 




9 9.5 10 10.5 11 

MR 



Figure 2: as a function of MR for MhR = 1-6 and s = 0. 

If s 7^ in this scenario the lowest eigenvalue is a Majorana mass given by 

~ Y^^v'^Mh [2ce-2'^l^*^l^ + coth(7rMi?) - ccoth(7rcMi?)] , (5.3) 

which in general will be too large, unless there is a strong suppression of the 5D Yukawa 
constants as y^,^ ~ 10~^, which is similar to the electron Yukawa coupling in this kind of 
models ije. 

Yet another possibility could be to localize the lowest mode of the RH neutrino towards 
y = 0, corresponding to cM > 0. Considering now MR ^ 1 the mass eigenvalue is proportional 
to (c + sign(M) + 2 e-^'^l^^l^ - 2c e"^'''^*^^). Then by choosing c = -sign(M) we could be left 
with an exponentially suppressed Majorana mass. However this value of c is not consistent with 
the initial hypothesis cM > 0. In fact the smallness of the Majorana eigenmass is naturally 
achieved with a different localization of the quark and lepton sector within the SM as we will 
see in the next section. 

-'^^As it is shown in Ref. [30] for values of MhR near MhR ~ 1-6, when supersymmetry is gfobaUy broken 
by Scherk-Schwarz boundary conditions the spectrum of the Higgs presents a tachyon at the tree level, which 
partially cancels the positive one-loop radiative correction to the Higgs mass due to the gauge coupling and 
allows EWSB to take place at the two-loop level with a modest amount of fine tuning. 
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5.2 Majorana mass 



The main obstruction to get a small Majorana mass eigenvalue out of the effective mass 
matrix (4.25) for the s 7^ case is that the Yukawa couplings of the higher KK modes are 
not suppressed if the SM matter is located on the same boundary where the Higgs localizes 
towards. However by allowing the Standard Model matter to be split into different branes 
for instance quarks localized in the same boundary (quark brane) where the Higgs localizes 
towards, and leptons localized in the opposite boundary (lepton brane) as Fig. 3 shows, then 
by means of the small exponential wave function factor e/^, the whole tower of effective Yukawa 
couplings will be exponentially suppressed by the Higgs localization and so /i will be. 




Figure 3: Bulk and brane matter distribution for a neutrino Majorana-like mass. The 
propagates in the bulk with mass M. 



Now the effective mass matrix is analogous to the previous case except for the global exponential 
suppression in the Dirac couplings, namely, fi ejj fi with en = g-'^^^w^. In addition we will 
assume the lowest mode of the RH neutrino to localize towards the leptonic brane, i.e. yj^ = 0, 

"'^^Such a splitting may find its justification within the context of intersecting branes in string theory. See 
e.g. [3f] and references therein. 
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corresponding thus to cM < 0. We then find that the lowest neutrino Majorana mass eigenvalue 
is given by 

m^J = ej^v^Y^Mn^-^ [2c + coth(7rMi?) - ccoth{c7iMR)] . (5.4) 

Notice the almost independence on the RH neutrino bulk mass M. Albeit its presence is 
absolutely necessary to provide the existence of a lowest Majorana mass eigenvalue it is 
shielded by the higher RH neutrino modes. In Fig. 4 we plot m^^ as a function of log^^g Vt^i 
fixed values of c and MR. 




2.6 2.8 3 3.23.43.63.8 4 

- logio y>yi 

Figure 4: Majorana neutrino mass, m^J, as a function of — log^g Vfi for c = —1/2, M^R = 1.6 and 
MR = h. 

From Fig. 4 we can see that generically < 1 eV implies ^ 10^^. 
The charged leptons, on the other hand, have masses 

mi ~ vyi^/MnReH , (5.5) 

By fixing in this scenario the Higgs localizing mass to its previous value MhR = 1.6 we can 
predict the correct value of the r mass [32] by means of a 5D Yukawa coupling ~ 1 while 
yi ~ mi/mr for the first two generations {i = e,fi). 

An interesting particular case arises here. Given that cM is negative, in the limit when 
\cMR\ > 1 Eq. (5.4) reads 

~ Y;^yelMH^-^ [3c + sign(M) + 2sign(M) e^^'^l^'l^ + 2c e^^^l^^^l^ . (5.6) 

^^In case of vanishing M we would be left with a lowest Dirac eigenvalue or, at most, with two almost 
degenerate Majorana eigenstates. 
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Considering now the particular value c = — |sign(M) we find a doubly suppressed Majorana 
mass eigenvalue. For instance for the case M > and c = —1/3 one gets 

< - ^yyele-h'^^ . (5.7) 

which has a doubly suppressed exponential behavior both from MhR and MR. In that case 
one can get tiny Majorana neutrino masses from the localization of the zero mode of vr for 
0{1) values of the 5D Yukawa couplings y^,^. This is shown in Fig. 5 where the Majorana mass 
m^^ is plotted versus \MR\ for = 1 and c = ±1/3. 




6 6.5 7 7.5 8 8.5 9 

\MR\ 



Figure 5: Neutrino Majorana mass, m^^, as a function of \MR\ for M^R = 1-6 and c = 1/3 (upper 
curve), c = — 1/3 (lower curve). 



6 Conclusions 

In this paper we have investigated the nature and size of the LH neutrino masses in a super- 
symmetric 5D model compactified in the space A^4 x / where / = [0, ttR] is a finite interval 
of length vri? and R the compactification radius. Quarks and leptons are localized on some of 
the boundaries and the gauge and Higgs sector propagate in the bulk of the fifth dimension. 
In this model, as we have found in previous works, supersymmetry can be globally broken 
by a Scherk-Schwarz twist giving a mass to gauginos (and gravitinos) which is transmitted 
by one-loop radiative corrections to squarks localized on the boundaries, and provides a very 
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characteristic pattern for the supersymmetric spectrum which could be easily identified exper- 
imentally at LHC whenever it is found [27]. Furthermore electroweak breaking proceeds by 
radiative corrections [30] providing at low energy an MSSM-like 4D model. 

We have found that the nature of the LH neutrino mass depends on the relative orientation 
of the boundary terms, given by the vector s, with respect to the bulk mass term characterized 
by the vector Mp. In the generic case of arbitrary orientations lepton number is violated and 
neutrinos are Majorana fermions. In the particular case where vectors s and p are aligned (or 
anti-aligned) there is an extra global U{1) symmetry which prevents lepton number breaking 
and LH neutrinos get a Dirac mass. 

As for the size of the LH neutrino masses, in the case of Dirac neutrinos, the smallness 
of neutrino masses is provided by the bulk mass M which should localize the zero mode of 
RH neutrinos towards the opposite boundary to that where SM fermions are localized: sub-eV 
masses are obtained for MR ~ 10. In the cases where neutrino masses are Majorana, in order 
to avoid introducing too small neutrino (dimensional) Yukawa couplings we have to localize 
quarks and leptons on different boundaries. In that case sub-eV neutrino masses are provided 
in general for Yukawa couplings ~ 10~^R. However for particular values of s ■ p there is a 
cancellation in the neutrino mass matrix in the limit of large localizing masses and the correct 
values ~ meV simply require MR ~ 10. 

Following the lines of our calculation it should be easy to describe textures of RH neutrino 
masses describing the different patterns for LH neutrino masses and mixings (see e.g. [33]). It 
should be enough to introduce the corresponding three-by-three mass matrices in the 5D bulk 
and boundaries and to carry on the parallel calculation. In the cases where Dirac or Majorana 
neutrino masses are controlled by the localizing masses of the RH neutrinos, since the former 
depend exponentially on the latter a modest change in the corresponding RH mass eigenvalues 
should be able to describe realistic neutrino spectra. Of course since the RH neutrino masses 
are an input in our theory, even if correct spectra do not require to fine-tune any parameters, 
we should not call this a "solution to the neutrino mass problem" until some more fundamental 
theory (e.g. string theory) would give us the correct values for the heavy masses. In fact this 
was not the aim of our work but rather a classification of the different solutions of the 5D 
theory providing realistic spectra for neutrino masses and yielding hints for possible future 
discoveries at LHC. 
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